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1 Introduction to Rotor-Router Model 

Informally a random walk is a process in which a particle takes a succes- 
sive random steps. A pseudo random walk appears to be a random walk, 
but it is not. It exhibits randomness in some sense but is obtained from a 
deterministic procedure. One of the popular random walk is the one that is 
executed on the lattice points. In the further paragraphs you will be able 
to notice that, the walks executed by the particles in Rotor- Router Model 
together can be viewed as a pseudo random walk. 

The need to develop the Rotor-Router Model arose during the study of 
Internal Diffusion-limited Aggregation (IDLA) and abelian sandpile models. 
IDLA is a model in which a particle executes a random walk in the lattice 
Z"^ beginning at the origin and ends at the point when it first reaches a point 
outside the set A which contains origin. Now the end point of this random 
walk is adjoined to the set A, and this process is iterated. The fact that after 
n walks have been executed, the resulting region A rescaled by a factor n^^'^ 
can be approximated by an Euclidian sphere in M*^ as n — y oo was proved 
by Lawler(t4j, 1992). Whereas in an abelian sandpile model, introduced by 
Bak et al.(|n], 1988) and studied by Dhar([8j, 1990) and Bonjour, et al.([7], 
1991), each of the points in A were assigned a non-negative integer, which 
represent the number of grains of sand at that point. Origin is assigned 
with a single grain of sand and at each step, every point that has at least 
2d grains, ejects a grain to each of its neighbours. The process is repeated 
until every point contains at most 2d — 1 grains. Many of the fundamen- 
tal mathematical properties which remain conjectural in IDLA can actually 
be solved in abelian sandpile model. Whereas no asymptotic results were 
known for sandpile models and it also appears likely that its asymptotics 
are not spherical. Though they have a common abelian property, there are 
considerable differences between IDLA and sandpile model. So a need arose 
to develop a model that is intermediate to both. 

Rotor- Router Model introduced by Jim Propp(pP) can be visualized as a 
deterministic analogue of IDLA. First a cyclic ordering of the 2 directions 
that are along the co-ordinate axes in Z*^ is fixed. A 'rotor' is placed at 
each point in A which specifies one of the 2d directions. When a particle 
is passed through a rotor it changes its direction to the direction that is 
following immediately after its previous direction in the cyclic ordering. A 
particle starts at origin and is routed through these points until it reaches a 
point that is not in the initial set A. Once the particle reaches to a point 



outside the set A, the same point is embedded to this set placing a rotor 
having an initial direction same as that of the origin. Then a new particle will 
be starting from the origin and follows the same procedure in new larger set. 
The process will be continued in this manner. 

These kind of models are very useful for computational purposes. They 
converge much faster than the conventional random processes. 

2 Rotor walk as a substitute to Random walk 

Consider the sequence {1,0,1,0,1,0,....}. This can be seen as a realiza- 
tion of the sequence of iid random variables {Xi} which takes values or 1 
with equal probability. Now consider the case of simple random walk in Z*^. 
Suppose if we have a collection of sequences of iid random variables at each 
lattice point, where each random variable assumes all the 2d directions with 
same probability. Now start at the origin and execute the random walk in 
the following manner. Whenever the particle reaches a point in Z'^ for the n*^ 
time, look at the n*'* random variable of the corresponding sequence present 
at that point. Based on the direction obtained from that random variable 
make the next move. In this way a simple random walk on Z'^ is executed. 

Now instead of the sequences {Xj}, one can use the sequences of the form 

{eo, ei, .., e2d-i, Cq, ei, .., 62^-1, Cq, 62, .., 62^-1, }, where {cq, ei, ...., 62^-1} are 

the natural 2d directions in Z*^. These sequences can be viewed as a realiza- 
tions for the sequences {Xj}. The walk that is obtained now is a rotor walk. 
In this way one can see the rotor walk on Z'^ as a pseudo random walk on 
Z*^. Rotor walk in general exhibits the properties of the random walk. 

3 Mathematical formulation of Rotor Walk 

Let Ed = {eo, ..., 62^-1}, be the set of 2d directions that are along the co- 
ordinate axes in Z"^ in that given order. A state of a Rotor-Router process 
is represented by a pair {x,l), where a; G Z'^ represents the location of the 
particle, and / : Z'^ — > {0, 1, ...2d — 1} indicates the direction of the rotor at 
each point. Define g{x, I) = {x + ei(^x), Ix), where l{x) is the labeling given by 

, , ,x f Hx') + Ifmod 2n) if x' = x: ,_,. 

Ux) = <! ^ ^ ,1. 1 

l{x ), otherwise. 



This explains the position to which the particle is moved and also the 
direction to which the rotor at that point is shifted. Composing the func- 
tion g to itself at the state (0, /) we can obtain the path of the particle as 
(0,xi,X2, ...) is a lattice path in l/, beginning at the origin. We denote this 
path by p = p{l)- Earlier we have assumed that the particle reaches outside 
the finite region containing the initial point with in finitely many steps. It 
will be proved in the following lemma. 

Lemma: Let A he the finite set containing the origin. Then the lattice 
path p{l) leaves the region A in finitely many steps. 

Proof: Suppose not, then 3 a lattice point in A such that the particle is 
visited infinitely many times. And for any point which is visited infinitely 
many times, the same true for its neighbouring lattice points. So any point 
that is connected with the point that is visited infinitely many times is also 
visited infinitely many times. Inducting the above argument along the path 
from origin to the point outside A, we obtain that the particle will reached 
outside the region A. But it contradicts our assumption. Hence proved. 

Thus the entire process in the Rotor-Router Walk can be interpreted math- 
ematically by the above terminology used by Levine([T]). 

Recurrent States: 

We say a state in !/■ is ^ recurrent, if a particle started at origin visits that 
state infinitely often. It is easy to see that if a state in Z*^ is recurrent then 
every other state in Z*^ is also recurrent. So the initial rotor configuration 
for which every state is recurrent is called 'Recurrent Rotor Configuration'. 
So if a state is not recurrent then we call it a 'Transient state',i.e. a particle 
starting at origin visits this state only finitely many times, so is every other 
state. Hence the initial rotor configuration which is not recurrent is a 'Tran- 
sient Rotor Configuration' 



An open problem from Holroyd and Propp([10]) 

We know that in a classical random walk, every state in Z'^ is recurrent 
for d < 2 and is transient for d > 3. We can easily find a several examples 
of initial rotor configurations for which it is recurrent in Z'^((i < 2), transient 
in Z°'((i < 2) and transient in Z'^{d > 3). In this context Holroyd and 
Propp([IO]) posed a question that for a Rotor walk on Z'^ with d > 3, does 
there exist a rotor configuration which is recurrent. 



4 Recurrent Rotor Configuration in Z^ 

Here in this section we would like to suggest an initial Rotor configuration 
so that it exhibits the recurrence property. We will propose a rotor config- 
uration in such a way that whenever the particle visits any site for the first 
time will be sent towards the origin. Mathematically it can be defined as 
follows. 

Let {xi, X2, ..., Xd) be any point in Z'^. Let cq, 62, ..., Cd-i be the unit vectors 
along the positive Xi, X2, ..., X^, axes respectively and e^, e^+i, ..., 62^-1 be the 
unit vectors along the negative Xi,X2, ■.■,Xd axes respectively. Now define 
the initial rotor configuration as follows: 

'' e2d-i if |a;i| > 'max{\x2\, \x3\, ..., 1x^1} and xi < 

ei if \x2\ > max{\xi\, \x3\, ..., \xd\} and a;2 < 

Crf-i if \xd\ > max{\xi\, \x2\, ■■■, \xd-i\} and a;^ < 

Cd if |a:i| > max{\x2\, {x^l, ..., 1x^1} and Xi > 

Cd+i if |x2| > max{|xi|, |x3|, ..., \xd\} and X2 > 

e2d-2 if \xd\ > max{|xi|, |x2|, ..., |xrf_i|} and x^ > 

ei in all other possible cases 

Also the sequence of the directions in which they rotate is (cq, ei, ..., e^-i, e^, e^+i, ..., e2d-i] 
The above will give us the initial rotor configuration in Z'^. 

We claim that the above described rotor configuration is recurrent. 

The case d=3 

For the case d = 3, we have a strong numerical evidence to show that the 
above claim we made is true. Define the box to be ]B[0,n] = [— n,n]'^. It 
is observed that the particle is visiting the origin 6n+l times before it first 
leaves the box B[0,n]. This has been checked for all n < 250. It is also 
observed that for any point after some time, the direction assigned to it is 
same whenever the particle exits the box ]B[— n,n]. The above mentioned 
facts clearly support our claim. 



/(xi,X2,...,Xd) = < 



5 What Next? 

If our claim gets proved, then one can attempt to classify the rotor con- 
figurations into transient and recurrent classes. Attempts can also be made 
to characterize these configuration by a simple rules. Studying these classes 
separately in detail may yield better results and can used as a replacement 
for classical probability models involving random walks. The transient con- 
figurations in l/'{d < 2) and recurrent configurations in U^^d > 3) are the 
analogues of measure-zero sets in classical random walk. Hence the study 
about these in detail can be used to model the conditional events of these 
measure-zero sets. 



6 Acknowledgements 

I thank Dr.Arni S.R. Srinivasa Rao for exposing me to this emerging field, 
and also for guiding me in this project for the past few months. I also thank 
all of my friends who helped me in this project. 



References 



[1] Levine, L.,"The Rotor-Router Model," | arXiv. ■math/04 094 0^ 1 
[math. CO]. 

[2] Levine, L. and Peres, Y., "Spherical Asymptotics for the Rotor- Router 
Models in Z'^," Indiana Umv. Math Journal 57 (2008), no. 1, 431-439. 

[3] Alexander E. Holroyd, Lionel Levine, Karola Meszaros, Yuval Peres, 
James Propp and David B. Wilson, "Chip-Firing and Rotor- Routing on 



Directed Graphs" arXiv:0801.3306v3 [math. CO]. 



[4] Lawler, G. F., M. Bramson and D. Griffeath, "Internal diffusion-limited 
aggregation," Annals of Probability 20 (1992), no. 4, 2117-2140. 

[5] Lawler, G. F., Intersections of Random Walks, Birkhauser, 1996. 

[6] Bak, P., T. Chao and K. Wiesenfeld. "Self-organized criticality," Phys. 
Rev. A 38 (1988), no. 1, 364-374. 

[7] Bjorner, A., L. Lovasz and P. Shor, "Chip-firing games on graphs," 
European J. Combinatorics 12 (1991), 283291. 



[8] Dhar, D. "Self-organized critical state of sandpile automaton models," 
Phys. Rev. Lett. 64 (1990), 16131616. 

[9] Sri Mallikarjun J., Srinivasa Rao S.R. Arni, "Generating Rotor- Router 
Models by a mapping on a Rotor- Router Model" M.Stat II year Project 
Report. 2008-09 

[10] Alexander E. Holroyd and James Propp, "Rotor Walks and Markov 
Chains" \arXiv:0904.450% j3 [math. PR]. 



